Bifurcations in Resonance Widths of an Open Bose-Hubbard Dimer 
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We investigate the structure of resonance widths of a Bose-Hubbard Dimer with inter-site hopping 
amplitude k, which is coupled to continuum at one of the sites with strength 7. Using an effective 
non-Hermitian Hamiltonian formalism, we show that by varying the on-site interaction term \ the 
resonances undergo consequent bifurcations. For A = fe/7 > 0.5, the bifurcation points follow a 
scaling law x-m = XmN/k = fhirn — 0.5/A), where N is the number of bosons. For the function /a 
two different A dependences are found around the minimum and the maximum bifurcation point. 
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The interplay of intrinsic dynamics with coupling to 
the continuum, radiation field or to any other external 
influence, such as measurement, is an important subject 
for various branches of modern physics: quantum chaos 
nanoscale devices 0, quantum optics 0], theory of 
decoherence and quantum computing Q up to nuclear 
[j| , atomic and molecular physics |(J are some of the ar- 
eas, seemingly far remote from each other, that boosted 
the research in open systems. These systems are often 
described within the "effective non-Hermitian Hamilto- 
nian" formalism 0. The eigenvalues of the effective 
Hamiltonian are complex £ n = E n — iT n /2, with the non- 
zero imaginary part, describing the rate with which an 
eigenstate of the open system (termed resonance state) 
decays in time due to the coupling to the " outside world" . 

With the advent of sophisticated techniques for trap- 
ping and transporting ultracold atoms, resonance states 
become an issue also in the context of Bose Einstein Con- 
densates (BEC). Examples include "atom lasers" [a, op- 
tical tweezers Q and magnetic "conveyor belts" g, as 
well as microtraps on "atom chips" \\(\ , which were sug- 
gested as potential building blocks for quantum informa- 
tion processing In this context, the aim is to under- 
stand the scattering and decay properties of a condensate 
leaking out of the trap via tunneling. 

First theoretical investigations of the waveguide scat- 
tering of a BEC through a double barrier potential, act- 
ing like a quantum dot for the atoms, reveal intriguing 
nonlinear effects associated with the internal resonances 
of the quantum dot ^2- Complementary studies on 
the decay of a condensate from one open trapping po- 
tential were reported in 0, Q, 0], while a first at- 
tempt to study the current relaxation of an ultracold 
BEC periodically driven with a standing wave of laser 
light was performed in [l^. The theoretical approaches 
used in these studies were based on the nonlinear Gross- 
Pitaevskii equation describing the mean-field dynamics 
of the condensate. 

However, a full quantum many-body study of decay- 
ing properties of a BEC is lacking. If the hopping am- 



plitude between nearby traps is small compared with the 
excitation energies in higher bands, one can use a Bose- 
Hubbard Hamiltonian (BHH) to describe the condensate 
quantum mechanically. This is the simplest non-trivial 
model that describes interacting bosons on a lattice, and 
incorporates the competition between kinetic and inter- 
action energy of the bosonic system. In this context, the 
closed two-site system (dimer) has been analyzed thor- 
oughly and many interesting phenomena have been found 
[H EH EES Ell- The richness of the results, provides 
a motivation to go beyond the closed dimer and consider 
new scenarios where even richer behavior should be ob- 
served. In this respect, the open quantum dimer promises 
new exciting opportunities, since the coupling to contin- 
uum leads to an interesting interplay with the internal 
dynamics. 

In this paper we address such a scenario and study the 
resonance linewidths T of a BHH coupled to continuum 
from one of the two sites. One possible realization of 
this setup is two coupled bosonic traps where tunneling 
to continuum is imposed at one site. Another motiva- 
tion comes from exciton transfer in molecular aggregates 
in which guest molecules (traps) are introduced intersti- 
tially and excitons are captured once they appear inside 
the sphere of influence of the traps [22l- Using an ef- 
fective Hamiltonian formalism, we show the appearance 
of consequent bifurcations in the resonance widths T n as 
we change the on-site interaction strength \. The reso- 
nances follow a scaling law 
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r = * A (x), x = X N/k (l) 



where N is the total number of interacting bosons, fc the 
hopping amplitude, 7 the coupling to continuum, and 
A = fc/7. In particular, we reveal that if the dimension- 
less ratio A = fc/7 is larger than A* = 0.5, then for N 
interacting bosons there are N/2 (2^| bifurcation points 
Xmi m = !)•••) A/2 which follow a scaling law 



Xr, 



,N 



/a (m) > M : 



m- 0.5/A 
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where /a represent a set of one parameter scaling func- 
tions which depend on A. Two different A— dependencies 
are found for the upper and for the lower part of the bi- 
furcation spectrum. The former is A-independent and 
thus universal while for the latter we have 



x-m/x* = g(v), 



x 



exp(-1.55/A) 



(3) 



where g(-) is a universal (A-indcpcndcnt) scaling func- 
tion. Our results i|H2l3|) are confirmed numerically and 
are supported by theoretical arguments. 

A usual way to treat the coupling of bound states with 
continuum is to exclude the continuum degrees of free- 
dom from consideration by the introduction of an ef- 
fective Hamiltonian H eS acting within the subspace of 
bound states and implicitly taking into account their in- 
teraction with the continuum p|. In the case of a dimer 
trap, it is natural to assume that the probability to es- 
cape in the continuum is proportional to the number of 
the bosons located at the site (sa y the second one) which 
is coupled to the continuum 24]. Following this line of 
argumentation we come out with an H given by 



£ eff = |X>(n;-l) 



; = 1 



i-fb 



(4) 



where \ = Anti a s V B { { /mo describes the interaction be- 
tween two atoms on a single site (V ott is the effective 
mode volume of each site, mo is the atomic mass, and 
a s is the s-wave scattering length of atoms which can be 
cither positive or negative and can be controlled exper- 
imentally). The operators = b\bi count the number 
of bosons at site i; the annihilation and creation opera- 
tors hi and b\ obey the canonical commutation relations 
\bi, b\] = Si , . We note that similar effective Hamiltonian 
was used to study a leaking dim er 1251 in the presence 
of pumping or driving light field |26j. The Hamiltonian 
Q in the fixed (instantaneous) basis of the eigenstates 
of TV = X)i=i reads 



H„ 



+ (TV - nf — TV] — i~f(N - n) 



n 



-ky/n(N + 1 - n) 
-k^(n+l)(N -n) 



:n = m + 
: n = m — 
: else. 

(5) 

where n,m = 0, 1, . . . , TV. We are interested in its (in- 
stantaneous) complex eigenvalues £ n = E n — i^f-, where 
E n and Y n are the position and the width of the reso- 
nances, respectively. The subindex n indicates the num- 
ber of bosons on the first (non-dissipative) site. Accord- 
ingly N — n bosons are on the second (dissipative) site. 
The rank of H eS is related to the number of particles as 
Af = N + 1. In our numerical analysis we consider traps 
with number of particles TV = 10 up to 500. 

To gain some insight in the parametric evolution of res- 
onance widths, let us first consider the two limiting cases: 




FIG. 1: Bifurcation diagram for the rescaled resonance 
widths r, vs. the interaction strength \ f° r three representa- 
tive coupling ratios: (a) A = 1, (b) A = 0.5 and (c) A = 0.4. 
Circles (o) and diamonds (o) correspond to k = 1 and 0.1 
respectively. In this particular example we consider TV = 10. 
The solid lines are drawn for the eye and show the flow of 
resonance widths as \ 1& changing. The stars (*) indicate 
the predicted asymptotic values of the resonance widths (see 
Eq. 100). The dashed (blue) line is the "classical" result 



vanishing interaction \ — and very strong interaction 
X> fc. In the former case one finds that the eigenvalues 
£ n of the effective Hamiltonian are given by 



1 



^)V A2 -i 



(6) 



From this expression one can see that for A > 1/2 the 
resonance widths are all degenerate and equal to TV7 
(see Fig. QJi,b). Instead, for A < 1/2 the eigenvalues 
£ n are imaginary and we do not have any degeneracies 
for x = (see Fig. QJ:). We have checked that whenever 
the widths show degeneracy the corresponding positions 
of resonances E n are not degenerate, exhibiting a "dual" 
behavior as far as bifurcation points are concerned. 
In the opposite limit of x ^ 1 one can neglect the 
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inter-site hopping term and thus, the Hamiltonian (0J 
becomes diagonal leading to 

£n = N,(^ [n* + (N-n?-N]-i(l-%)). 

. (7) 

Therefore, we obtain N equidistant resonance widths 
r„ = 2j(N — n) corresponding to n particles being 
trapped on the first site (which resembles the resonance 
trapping phenomena known for the non-interacting parti- 
cles |27| ) . Naively one could think that the initial degen- 
eracy of the resonance widths for A > 1/2 can be lifted 
completely by arbitrary small interaction x and the reso- 
nance widths flow continuously to their limiting values by 
increasing x- Below we will see that surprisingly enough 
this is not the case and the degeneracy is reduced each 
time by 2 at N/2 discrete points. 

In Fig. ^ we report the resonance widths for three rep- 
resentative values of the coupling ratio A = 1, 0.5 and 0.4, 
as a function of the interaction strength. Our numerical 
data are rescaled according to Eq. The points corre- 
sponding to the same value of A (but different values of 
k) fall onto the same smooth curve. In the same figure 
we report the asymptotic values given by Eqs. 10 0). 

In order to shed some light on the emergence of the 
bifurcations and the scaling ansatz Q we adopt a semi- 
classical point of view for Hamiltonian justified in the 
limit of N 3> 1. The corresponding classical eigenvalue 
problem reads 

Eai = x I ai I 2 ai - ka 2 , (8) 
£a 2 = x I a 2 | 2 02 - ka\ - i^a 2 , 

which can be solved exactly. We find that for x < Xcr = 
V4 — A~ 2 there are two "trivial" solutions with constant 
(%- independent) value of T = — > T = 1 corresponding 
to the particles distributed equally between the two sites 
| ai | 2 = \a 2 \ 2 = N/2. For x > Xcr two new non-trivial 
solutions appear with =jN(l± \/l - 4/(x 2 + A" 2 )) 
corresponding to the non-equal occupations \af\ 2 = 
Y ( 1 T a/1 - 4/(x 2 + A -2 )) . The classical results are pre- 
sented also in Fig. ^ ( see dashed lines). This pitchfork 
bifurcation is the classical analog of the quantum bifur- 
cations discussed above. 

We see further that for A > 0.5 resonances are initially 
degenerate, as it is predicted by Eq.JfJJ). As x increases 
this degeneracy is lifted by consequent bifurcations. Con- 
trary, for A < 0.5 we do not observe any initial degen- 
eracy. We note that in the regime A < 0.5 the instan- 
taneous approximation is questionable since the hopping 
time Tfc ~ 1/k is large in comparison with the character- 
istic decay time r 7 ~ 1/7. Therefore, in the rest of the 
paper we will concentrate in the opposite case A > 0.5 
and analyze the sequences of bifurcation points Xm which 
we have extracted from our numerical simulations. They 
were identified as the points after which two degenerate 
resonance widths T differ by more than C/N, where we 




FIG. 2: (Color online) The bifurcation points Xm plotted 
against the scaled index \i for representative values of the 
dimensionless ratio A = 0.65, 1,10 6 . The symbols corre- 
spond to different boson numbers iV and different hopping 
amplitudes k: N = 10, k = 1 (black), N = 30, k = 10 
(red), N = 50, k = 20 (blue). The symbol A corresponds 
to TV = 20, k — 0.1. The solid line indicates the large N limit 
of the case A = 10 6 , while the dashed line corresponds to the 
value A = 0.5 and a system size of N = 50. In the inset we 
show a magnification of the main figure. 

used C = 0.01. In Fig. [5] we report our findings by mak- 
ing use of the scaling proposed in Eq. J5J). An excellent 
agreement is evident. We have checked that our results 
remain qualitatively the same if we apply other similar 
criteria. A side remark is that the resonance positions 
E n (not shown) bifurcate at the same points Xm- 

Now we focus on the dependence of the scaling func- 
tion /a on the parameter A. Our study has revealed the 
existence of a sort of transition separating two regimes 
characterized by different scaling properties. Namely, we 
found that at large /^-values (i.e. large x) the various A- 
curves fall nicely one on top of another, even for A = A*, 
as clearly demonstrated in the inset of Fig. [3 From our 
data we were able to estimate that this scaling of the 
upper part of the spectrum holds for fi > 0.35. Notice 
that in the limit of N — > 00 at the fixed ratio m/N the 
scaling parameter fi ps m/N i.e. is independent of A. 
In other words the bifurcations with large m values take 
place between states that are not affected by the inter- 
site coupling k and the coupling to continuum 7. Indeed 
these states correspond to almost equidistributed number 
of bosons between the two traps. 

On the other hand, for /i < 0.3 we see that all curves, 
corresponding to different values of A have the same func- 
tional form, albeit being shifted downwards with respect 
to each other (Fig.[21). The curves x do coincide, however, 
when rescaling the bifurcation spectra as Xm/x* where 
X* is a scaling parameter (Fig. |3 main part). We find 
that the scaling parameter x* depends only on A (Fig. |3 
inset) resulting in a universal one-parameter scaling for 
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FIG. 3: (Color online) The rescaled bifurcation points Xm/x* 
for three values of A = 0.65, 1, 10 6 . The symbols correspond 
to different boson numbers iV and different hopping ampli- 
tudes k (see Fig. 2) . The solid line indicates the large N limit 
of the case A = 10 6 . Inset: the extracted scaling parameter 
X* as a function of the dimensionless ratio A . 



the lower part of bifurcation points according to Eq. J3J). 
The best fitting indicates that x*(A) ps exp(— 1.55A). We 
note that the scaling region becomes increasingly small 
as A — > A* and breaks down totally at A = A*. This 
is the limit where the resonance widths do not show any 
degeneracy for x = 0. 

In conclusion we have studied the dependence of res- 
onance widths of an open dimer BHH system on the 
strength of the on-site interaction and found a rich bifur- 
cation behavior. Our results have immediate implications 
in the dynamics of an open dimer. Specifically, we ex- 
pect that the decay of the bosons is not given by a simple 
exponential as suggested by the standard rate equation 
TV = — FN with r ~ 7 being a constant. Instead, our 
analysis indicate a more complicated decay since now F 
depends on the remaining particles inside the trap (see 
Eq.JU and Fig. 1) and thus one has to solve a non-linear 
rate equation. The analysis of dynamics of an open BHH 
system is beyond the scope of the current publication and 
will be the subject of a separate paper. 
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